Energy Spectra of Few-Electron Quantum Dots 

E Anisimovasf § and A MatulisJ || 

f Vilnius University, Sauletekio 9, 2054 Vilnius, Lithuania 

I Semiconductor Physics Institute, Gostauto 11, 2600 Vilnius, Lithuania 

Abstract. We present the renormalized perturbation series for the energy spectrum 
of the parabolic quantum dot with 2-5 electrons considering ground and the lowest 
excited states. The proper classification of asymptotic energy levels is performed and 
behaviour of energy levels from quantum to semiclassical regime is traced. Comparison 
between the present results and those of exact numerical Hamiltonian diagonalization 
shows a fair accuracy of the proposed method over the whole range of the electron- 
electron coupling constant and magnetic field values. The obtained results indicate 
that increasing of the number of electrons in a dot leads to more classic behaviour of 
the system. 



PACS numbers: 71.10.+x: 73.20.Dx: 03.65.-w 



Short title: Energy Spectra of Few- Electron Quantum Dots 
February 1, 2008 



§ Present address: Department of Theoretical Physics, University of Lund, Solvegatan 14A, S-223 62 
Lund, Sweden. 

II Electronic mail: matulis@uj.pfi.lt 



2 

1. Introduction 



Quantum dots have been studied extensively as low- dimension nanostructures 
interesting from the intrinsic theoretical point of view and due to their importance to 
modern electronic devices The most spectroscopy experiments were carried out on 
quantum dots with soft parabolic confinement potential 0. Those parabolic quantum 
dots attract the attention as simple model systems for electron-electron interaction 
studies. Although the parabolic quantum dots have some disadvantage due to Kohn 
theorem |^ which shows that only the center of mass motion usually can be excited (in 
absorption and even in tunneling), the confining potential in a dot very often is rather 
close to the parabolic one. That is why the parabolic dots can serve as a good starting 
point or the reference for the investigations of more sophisticated quantum dots. The 
small nonparabolicity usually helps to reveal peculiarities of the relative electron motion 
spetrum which is sensitive to the electron-electron interaction. 

Only in the case of quantum dot with two electrons the exact analytical solution of 
the eigenvalue problem for the parabolic dot can be obtained 0]. The most of theoretical 
investigations were mainly based either on some mathematical tricks (which is possible 
in the case of parabohc dot with 2 and 3 electrons P, |^), or on the straightforward 
diagonalization of the Coulomb interaction [|^. Although the latter gives the exact 
solution of the problem, due to time consuming calculations it can be applied only in the 
case of small number of electrons in a dot. That is why the construction of approximate 
solutions is of interest. One of them was demonstrated in |Q where in the case of high 
magnetic field rather accurate results were obtained using the moving Eckardt frame. 
Recently it was shown that some approximate technique based on the renormalized 
perturbation series can be successfully used for the energy spectrum estimations in a 
wide coupling constant range. The other technique based on Pade approximants was 



shown to be very successful in the case of parabolic dot with small number of electrons. 
Both of those approximations are based on some interpolating technique between the 
perturbation series in coupling constant and the asymptotic expansion valid in large 
coupling constant range. 

In the present paper in order to obtain the total picture of energy spectrum, 
including the ground and some excited states, we performed the proper classification of 
the asymptotic energy levels paying attention to the symmetry properties of the spin 
wave function part and constructed the renormalized perturbation series for the energy 
of quantum dots with 2-5 electrons based on rahter simple two term expansion in both 
limiting regions of the coupling constant. 

The present paper is organized as follows. In section 2 the problem is formulated and 
the perturbation series technique is outlined. The asymptotic series are formulated in 
section 3 and the spin wave function part for asymptotic region is constructed in section 
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4. The interpolating technique is described in section 5 and in the last section 6 the 
results for quantum dots with 2-5 electrons are presented and discussed. In appendix 
A the algebraic expression for Coulomb matrix element is given and in appendix B all 
coefficients which are necesssary for construction of the renormalized perturbation series 
are collected. 



2. Perturbation series 

Generally the system of N electrons in a quasi-2D quantum dot is described by the 
following Hamiltonian 



e 



fc2 2 



i=l 

Here m stands for the effective electron mass, e is the static dielectric constant, and 
ujq is the characteristic frequency of parabolic confinement potential. For the sake of 
convenience we shall scale the coordinates r aov (oq = ^/H/muJo) and measure 
the energies in dimensionless units of Huq. After this transformation we obtain the 
dimensionless Hamiltonian 

i=l i<j ' * ' 

with the single dimensionless coupling constant \ = a^/ as {as = tf? jrrK? is the effective 
Bohr radius) characterizing the energy spectrum of electrons in the parabolic quantum 
dot. 

The simplest way to solve the eigenvalue problem is just to use the perturbation 
series in coupling constant A 

E{X) = Eo + E^X. (3) 

The first coefficient is the eigenvalue of the unperturbed Hamiltonian Hq and can be 
presented as the sum of noninteracting electron energies 

N 

Eo = ^e{i), e{i) = 1 + \mi\ + 2ni (4) 
1=1 

where i = {ni,mi,Si) labels the one electron quantum number set. 

The first order correction Ei can be estimated in a standard way (see, for instance 



Tl| ) and presented as a sum of two-electron Coulomb matrix elements 

1 



{zj\V\j'z') = 6s,,sjs„s, I rfVi / ^ r(z|ri)0-(j>2)0(/|r2)0(zVi) (5) 
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where the one electron wave functions can be expressed in associated Laguerre 
polynomials 



(n, m\r, ip) 



2n\ 



27r(|m| + n)\ 



(6) 



For numerical estimation of the Coulomb matrix element (|^) we used the analogue of 
Girvin and Jach expression |T^ which is given in [Appendix A| . 

When calculating the first order correction the main problem is the degeneracy 
of zero order energy levels because the good quantum numbers of the total angular 
momentum M and the total spin S are not sufficient for classification of energy levels. 
That is why we used numerical diagonalization of the degenerate blocks of the above 
first order correction matrix. The obtained results - the coefficients Eq and Ei are 



collected in Appendix B 



3. Asymptotic expansion 

In the A ^ oo case the kinetic energy of the electron system is small if compared to the 
potential one. It is known that such systems are strongly correlated and show tendency 
for Wigner crystallization |0 . Monte Carlo results for classical particles in the parabolic 



dot ||I4[ show that the particles crystallize in some structure of concentric rings. When 
the number of electrons in a dot is small N < 5 they fall into a single ring on which 
the electrons are located equidistantly and separated by the angle a = 2ti/N. The 
quantum dynamics of such a Wigner crystal can be described by means of the rotation 
and vibration modes. The general theory of vibration modes is given in It is based 
on introducing the moving Eckardt frame and the diagonalization of the vibrational 
Hamiltonian in the harmonic approach. In the case of small number of electrons located 
on a single ring the vibrational Hamiltonian and the problem of its diagonalization can 
be essentially simplified by representing the electron coordinates as complex variables 
z = X + iy and turning them to the local electron coordinates on the ring by means the 
following transformation 

^„^(i^ + ^„)e^"^ (7) 

Here n labels the electrons on the ring. We shall follow the simpler way restricting our 
consideration by the following two term asymptotic expansion 

E{\) = CoA2/3 + ci. (8) 

The first term of it can be easily obtained minimizing the potential in (0). The 
minimization leads to the equilibrium ring radius R = (XAY^'^ (where A = 
^^Ji^ |4sin(Q;n/2)|^^) which corresponds to the main term CqA^^^ = l.^NR^ in the 
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asymptotic energy expansion. The next term in expansion (j^) corresponds to the ring 
oscillations. To obtain it we insert (0) into Hamiltonian and expand the potential into 
2;-series. It leads to the following second order correction to the equilibrium potential 

1 ^ N 



2 ^ ' ' 32A 

n=l n=l 



where 



[l _ gia{n-m)j ^ _ gm(m-n)j 



Q^an_^iam 4 siu^ [« (n - m) /2] 

Now let us apply the Fourier transform 

l=J]e^^-.. (9) 



where k = 0,1,. ..,A^ — 1 denotes the vibration modes. Inserting it into the above 
potential expression and subsequently into (0) after simple but laborious transformation 
we obtain the following expression for the asymptotic Hamiltonian with the second order 
potential corrections 

H = J2Hk = Yl {-^d^ + ^'^^'^^^ + M^kZ-k + zlz*A (10) 



where 



1 1 

^k = - + —Tj2\ sin(an/2)|-3 sin2[(A; + l)an/2], 



2 16A 

n=l 

3 



Bk = — — I sin(an/2)|-^[sin^(an/2) - sin^(A;an/2)]. 

n=l 

The Hamiltonian is composed of parts which describe the behaviour of both 
polarization of ±k modes as a system of four coupled harmonic oscillators. They can be 
diagonalized in a standard way just introducing the proper variables. Let us separate 
the real and imaginary parts of the initial variables, and introduce the new variables in 
the following way 

Zk = x + iy, z-k = u + iv. (11) 

In these variables every Hamiltonian (|T0|) part Hk splits into two terms. The first term 
is 

1 / 0'^ 0'^ \ 

H^^^ = -2 (9^ + + ^'^^^ + + + B-k)xu (12) 

while the other part HI can be obtained just replacing x y and u ~v. 
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Both Hf' have the same eigenvalues that can be found by the diagonahzation of 
the symmetric matrix 

which defines the quadratic form of the potential term, and leads to the following 
expression for the eigenfrequencies 

= \jAk + ± ^[Ak - A_k? + 4[fi, + (13) 

Thus, since the transformed Hamiltonian corresponds to the set of independent 
harmonic oscillators the corresponding vibration wave function $ can be constructed in 
a standard way introducing the "phonon" creation ^ and annihilation a^^ operators 

'^ = <A---«UlO). (14) 

The "vacuum" state |0) corresponds to the system ground state, and the symbol [3 = 1,2 
stands for polarization. 

For the interpolation purposes it is important to know the wave function symmetry 
properties with respect to the rotation of the electronic ring as a whole. Let us denote 
this operation by operator V. It is evident that the ground state is invariant under 
that operation, i. e. "PlO) = |0). As the above rotation is equivalent to the substitution 
n — i> n + 1 it follows from expression (|^) that Vzj. = exp{ik)zk. Due to the mixing 
of both ±k modes the operators ^ and wave function ( |T^ obey more complicated 
symmetry condition. However, because of the degeneracy of obtained vibration states 
it can be essentially simplified using the proper superposition of the eigenstates of both 
Hamiltonian parts Hj^ ' . So, we shall assume that the wave function superposition 
of this sort is constructed, and consequently the phonon operators obey the symmetry 
condition 

-Pala = (15) 

We shall use the above symmetry condition for the classification of obtained 
asymptotic vibration modes whose frequencies ulf^ are collected in table |l|. Two modes 
with eigenfrequencies uj[^'^ = 1 correspond to the center of mass motion. We do not take 
them into account as we consider the electron relative motion energy spectrum only. The 
mode k = corresponds to the 'breathing' mode for (3 = 1 and to the rotation of the 
ring as a whole for (3 = 2. The last case should be considered separately. Actually, in the 
harmonic approximation used the rotation mode is separated from the vibration modes. 
Its eigenfunction is exp(zM(y9) {ip is an angle determining rotation of the ring as a whole) 
and the corresponding energy is M"^ /2NR^ ~ A^^/^. It is seen that in the case of large 
angular momenta M this term becomes essential and thus it should be included into the 
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Table 1. The ring vibration frequencies and the asymptotic expansion cocfBcicnts 
N 2 3 4 5 



k /? 1 2 1 2 1 2 1 2 

1.7321 0.0 1.7321 0.0 1.7321 0.0 1.7321 0.0 

1 1.0 1.0 1.2247 1.0 1.3186 1.0 1.3719 1.0 

2 1.2247 1.0 1.4888 0.8852 1.6890 0.7274 

3 1.3186 1.0 1.6890 0.7274 

4 1.3719 1.0 

Co 1.1905 3.1201 5.8272 9.2801 

c\°^ 1.8660 3.0908 4.3716 5.6544 



main asymptotic energy term, and consequently it leads to the equilibrium ring radius 
dependence on the total angular momentum M (see ||^ for instance). However, the 
above term is of order A"^/^ and lies beyond the (|]) approximation. Hence considering 
the states with small M only we neglect this rotation mode correction (and the ring 
radius dependence on M as well) together with the anharmonicity corrections which are 
of the same order of magnitude. We should keep in mind, however, that there are many 
states labeled by different M values corresponding to every asymptotic vibration term. 
Thus, the first order correction is given by 

^ _ ^(0) , „ , ,{/3) 

where ci' is the contribution of the ground vibration state energy and the symbol 
rifc^^ stands for the filling factor of the corresponding phonon mode. Those coefficients 
together with first term coefficients cq are given in the last two lines of table |1|. 

Separating the rotation mode we shall represent the total wave function in the 
asymptotic region as a product of rotation, vibration and spin parts 

^ = e**^^$T. (16) 

According to the Pauli's principle the total wave function should be antisymmetric. 
In the asymptotic region the classical electron equilibrium positions are separated by 
high potential barriers. That is why we shall restrict ourselves only to the physically 
important permutation — the rotation V . As the rotation corresponds to the odd 
permutation in the case of the even number of electrons (and vice versa) the total wave 
function obeys the following symmetry condition 



p$ = (-l)^-i^. 



(17) 
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According to (^^-[T6|) the rotation and vibration parts obey the condition 



(18) 



The prime in summation indicates that the mode corresponding to the ring rotation is 
excluded. The different factors in the two above conditions should be compensated by 
the proper choice of the spin function part. 

4. Spin wave function part 

Following [§] we diagonalize the V operator in the 5*^ (total spin 2;-projection operator) 
eigenfunction space and construct the spin function part obeying the symmetry 
condition 



In the case of odd number N of electrons in a dot the symbol p stands for an arbitrary 
integer while in the case of the even it should be replaced hj p + 1/2. The integer 
parameter r which characterizes the spin function part symmetry properties should be 
defined for every particular case. In the case when the above procedure leads to the 
degeneracy of spin states the total spin operator S"^ should be diagonahzed additionally. 

The case of 2 electrons is a trivial one. There is the triplet function | |"f) with 
Sz = I corresponding to r = and the singlet function {| ti) — I iT)}/v^ with Sz = 
corresponding to r = 1. We shall not consider the functions with < as they give 
no extra information. 

The case of 3 electrons is considered in detail in [§]. There is quartet state 
corresponding to r = and two doublet states with r = ±1. 

Applying the same treatment to the system with 4 and 5 electrons we obtain the 
spin function parts with the symmetry properties of ([19|) type. All the necessary r 
values are collected in table ||. 

The above selection rule gives us the possible terms for every vibration state 
characterized by 7 value. The terms for the ground vibration state (7 = 0) are collected 
in the last column of table For the sake of brevity we use the notation analogous 
to that adopted in atomic spectroscopy. Thus the capital letters S, P, D, ■ ■ ■ indicate 
the corresponding M = 0, 1, 2, ■ ■ ■ values, and the left superscript symbol stands for 
multiplicity (25*+ 1). The terms corresponding to the excited vibration states with 7 7^ 
can be easily obtained from those presented in table ^ applying the shift M ^ M + 7. 
Further we shall label those excited terms with the additional right underscript symbol 
k indicating the modes of the phonons present in that state. 



(19) 



which together with (|T7p, (|TBp leads to the following selection rule 
M + T - 'J = Np. 



(20) 
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Table 2. r values for the various multiplets and the ground state terms. 
N S = 5=1/2 5=1 5 = 3/2 5 = 2 5 = 5/2 Terms 

2 1 - - - - ^S,^P,^D,^F, ■■■ 

3 - 1,2 - - - ^S,^P,^D,^F, ■■■ 

4 0,2 - 1,2,3 - - i'35, i-^^L), 3^, • • • 

5 - 0,1,2,3,4 - 1,2,3,4 - 2,65 2.4p 2,4^, 2,4^ 2,4^ .. . 



5. Interpolation 

In constructing tlie interpolating expression we shall follow where the perturbation 
series renormalization procedure is given in details. The main idea is replacing of (Q) 
by the generalized Hamiltonian 

ff^iEKM+^a+A^:-^, (21) 

i=l i<j ' * ^' 

and using the scaling relation for the generalized eigenvalue 

E(e,A) = ei?(i,Ar'/'). (22) 

The eigenvalue of the basic problem with Hamiltonian (^ is given as E{X) = E{1, A). 
The scaling relation defines the trajectory A = Aq^^'^^ in the ^A-plane as it is shown 



in figure Al. That trajectory enables to map the vertical dashed line ^ = 1 (where we 
need to obtain the eigenvalue of the basic problem) onto tilted line KP which we define 

as 

X = K(3, ir=(l+tan</))/2 

e = P(l-/3), P = (l + cot0)/2. ^ ' 

For instance, the point Aq is mapped to the point (^, A). 

Contrary to line ^ = 1 the line KP is located in the finite region of the ^A-plane 
(0 < /3 < 1), and that is why the energy E{^{P), X{P)) can be successfully expanded 
into /3-series. Adjusting the above series to the A-series (^ in the region A — > and to 
the asymptotic expansion (P) in the region A — > 00 we obtain the final expression for 
the renormalized series 

E={bo + hP + b2P' + hP')/{l-P) 

A = (2A) 

The expansion coefficients can be expressed in terms of those in (^ and (|]) as 

bo = Eq 
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61 = KP-^/^Ei - Eo 

62 = 7/3 K^/^p-^co + C1-E0- 2KP-'^''^Ei 

h= -A/?,K'^''^p-^co-ci + Eo + KP-'^/'^Ei. (25) 

As we restricted ourselves with the two term approximation in both (|^) and (§) 
expansions we introduced the additional adjustable parameter (as proposed in |jl5|) 
for improving the interpolating accuracy. The geometric meaning of that parameter is 
clear from figure — it indicates the rotation angle of the mapping line KP around 
the point M to which the energy E{\/2) is always mapped. The rotation changes the 
relative contributions of the A-series and the asymptotic expansion to the renormalized 
series. The value = corresponds to the full control of the renormalization procedure 
by the A-series and the value = 7r/2 to that of the asymptotic expansion. 

It was shown in [|T5| that usually the energy exhibits some plateau as a function of 
the parameter controlling the mapping line (0 in our case). When the number of terms 
in the A-series is increased that plateau shows the tendency to increase, too, thus giving 
a good reason for the controlling parameter choice. In figure ^2| the energy at the fixed 
point A = \/2 as a function of the parameter is shown by the solid line in the case 
of the first excited triplet state of the two electron system. As our A-series are rather 
short we have no prominent plateau. However, some curve twist is present, and we hope 
that it is a rather good idea to choose the controlling parameter somewhere between 
the Max and Min points. The difference of energy values corresponding to those points 
can serve as an estimation of the interpolation accuracy. We see that the accuracy is of 
about 1%. 

We used the inflection point of E{\/2) versus curve for the controlling parameter 
choice. The energy value obtained with that choice is shown by dashed line in figure [A^ . 
Comparing it with the exact result which is known for two electron case (dotted line) 
we see that the accuracy is even better than 1%. The same energy versus behaviour 
and the same accuracy was obtained for all other considered terms. The inflection point 
values of for all the considered terms are given in [Appendix B| in tables [Al| and 



They together with coefficients i?o, Ei and expressions (p^25) define the renormalized 
perturbation series and enable to obtain the energies in the whole range of A values. 

We should note that according to the same parameter values can be used for 
constructing the renormalized series in the case when the additional magnetic field B is 
applied in the direction perpendicular to the quantum dot plane. In that case we have 



E{\ B) = Vl + (cj,/2cuo)2E(A') + Muj,/2iUo + g*fXBBS,/nwo (26) 

and 



(27) 
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where the symbol = eB/mc stands for cyclotron frequency. Note that there is a 
misprint in (5.4) of [Q. The last term added in (p6D corresponds to the Zeeman energy 
where g* is the effective g factor and /i_B is the Bohr magneton. 



6. Results and discussion 

As an illustration of the introduced renormalization technique we present in figures ^ 



A6| the excitation energy versus A plot for the relative electron motion for the quantum 
dot with 2-5 electrons. The inclusion of the center of mass motion energy is a trivial 
one. One should just shift up the presented spectra by the some integer value, as the 
center of mass motion is not affected by the electron interaction. 

The main problem is to decide how the terms in both regions (A ^ and A —>■ oo) 
should be interconnected. Taking into account the fact known in quantum mechanics 
that two terms with the same good quantum numbers never cross each other we used 
the simplest term connection scheme. Starting from the ground state we connected 
successively term in the region A ^ to the lowest possible unconnected term with the 
same M and S quantum numbers in the asymptotic region. 

In figures [A3| and the energy plot for the system of two and three electrons are 



presented. The terms with the zero energy exceeding the ground state energy by 3 and 
2, correspondingly, are depicted. In the case of four and five electrons the energy plots 
are given in figures and The degeneracy of the states grows rapidly with the 
number of electrons and thus here we present only the two lowest terms tending to each 
of asymptotic vibration states. The asymptotic terms in the above figures are labeled by 
the numbers in parenthesis along right axis which indicate the k values of the phonons 
present in the state. The phonon state with polarization corresponding to the larger 
frequency is indicated by prime. 

It can clearly be seen how the term structure transforms from that of the 
weakly interacting quantum system, where the energy levels are governed by filling 
of noninteracting single electron states to the semiclassical strongly correlated system, 
where the energy levels correspond to the electronic ring vibration frequencies. Such 
switching of the dot behaviour can be explained by simple dimensional considerations. 
Namely, the three terms in Hamiltonian (|l|) can be estimated as h'^N/mR'^ , uujJqB^N and 
e^N{N — 1)/ eR. Equating them we define the critical A value A* ~ 1/ (A^ — 1) where the 
above mentioned switching occurs. This is in agreement with our energy dependencies. 
As the dimensionless coupling constant A is proportional to oq ~ the above 

estimation leads to some condition for confinement potential frequency showing whether 
the dot behaves like quantum system or the classical one. When the confinement 
potential is mainly caused by some positive neutralizing charge located on the rather 
distant gate the increase in confinement potential frequency ujq does not exceed A^ M^. 
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In that case increasing of the number of electrons in a dot leads to more classical system. 

As an illustration of the presented renormalized series technique in the case with 
magnetic field (according to (^)) we calculated the ground state energies vs. magnetic 
field for 2-5 electrons which are displayed in figures (a)-(d). The curve (a) for 
two electron system is in good agreement with that presented in ||^ (we choose the 
parameters as in 0] and also include the fictitious gate charge but do not subtract 
the zero point motion of the center of mass). It is seen from figures (b)-(d) 
that increasing the number of electrons the structure of the dependencies becomes less 
resolved and the peaks lose their sharpness. Thus we miss a tiny additional peak in 
the position marked by an arrow in figure [A7| (b) which is present in data of exact 
diagonalization. Comparing our results with those obtained in [§], say, for = 3, 
B = 20 T, S = 3/2, M = 3 we obtained ground energy value 77.5 meV and find a good 
coincidence which proves that the approximate separation of rotation and vibration 
modes works fairly well for small angular momenta. 

Concluding we would like to mention that the parabolic quantum dot due to the 
softness of the confinement potential is rather favourable object for constructing the 
interpolating expressions in the contrary to the hard wall dot Ijl^ where the accuracy 
is much worse. 
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Appendix A. Coulomb matrix element 



In |jT2| the algebraic expression for calculation of the two-particle matrix element of 
Coulomb interaction of electrons in the lowest Landau level (corresponding to rij = 
in our case) was used (see the comments in [0). In more general (n^ ^ 0) case the 
analogous Girvin and Jach expression can be derived employing the same center of 
mass and relative motion complex coordinate technique after expanding the associated 
Laguerre polynomials in sums and representing the terms in form 

In - 

here the symbols 7^ stand for some numbers 

7i = Ji + J4 + (|"^i| + mi)/2 + (|m4| - m4)/2, 

74 = Ji +J4 + (l^il -mi)/2 + (|m4| +m4)/2. (Al) 

Symbols ji are integer summation indices running form to rij. 72 and 73 can be 
obtained from (Al) replacing the indices 1 — > 2 and 4^3. The final expression reads 

1/2 

„ 1 

\12|V^|34) 6sj^^s4^S2,S3^mi+m2,m3+m4 
' 4 

n 



X 



i=l 



where 



ni n2 n-j 




(A2) 



E -EEEE 

(4)j=o ii=Oi2=Oi3=oi4=o 



G = E 



A 




Tn\(n — m)! 

Although the above expression involves terms with alternating signs and leads to 



numerical difficulties in the case of large quantum number values [0 , we found it rather 
convenient and accurate in the case of quantum dots with small number of electrons. 



Appendix B. Expansion coefficients 



Here the perturbation series (0) coefficients Eq and Ei, and the inflection point angle 
are presented for all terms under consideration. Those values together with the 
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asymptotic expansion (j|) given in table ^ and expressions (p^25) are sufficient for 
constructing the renormalized perturbation series. 

Table Al. Coefficients for quantum dot with N — 2 and = 3 electrons. 



N — 2 








/V — 3 








Term 


Eo 


El 




Term 


Eo 


El 






2 


1.2533 


0.3571 


2p 


4 


2.8200 


0.5871 


3p 


3 


0.6267 


1.2211 




5 


1.8800 


1.0722 




4 


0.4700 


1.3993 




5 


2.2325 


0.9358 


'So 


4 


0.9400 


0.7204 


'Si 


5 


2.5850 


0.5836 


3f 


5 


0.3917 


1.4651 




6 


1.7037 


1.1984 


3 p 


5 


0.5483 


1.3065 


'Pi 


6 


1.8212 


1.0783 










'Pi 


6 


1.8413 


1.0698 










^Fi 


6 


2.2618 


0.8991 










'Po 


6 


2.5942 


0.7067 
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Table A2. Coefficients for quantum dot with TV = 4 and N = 5 electrons. 











N — ^ 
1 V — 








Term 


Eo 


El 




Term 


Eo 


El 




35 


6 


5.0133 


0.6105 




8 


8 3032 


0.4913 




6 


5.2483 


0.5800 




9 


7.4416 


0.6649 




6 


5.4832 


0.5527 




9 


7.4812 


0.6607 


3p 


7 


4.3608 


8333 


^2 


9 


7.6765 


0.6409 


3F 


7 


4.6412 


0.7828 


'^So 

^2 


9 


7.7647 


6325 


1 


7 


4.6999 


0.6087 




9 


7.8234 


0.6270 


■*■ 2 


7 


4.8628 


0.6482 




9 


7.8284 


6265 




8 


3.5641 


1.0635 


^Do 

^2 


9 


8.1710 


5969 




8 


3.6816 


9920 


'^So 
^2 


9 


8.3125 


5263 




8 


3.9491 


8956 


4p 


10 


6.7149 


8303 




8 


3 9822 


0.9218 




10 


6 8596 


0.8110 


'Do 

^2 


8 


4.0215 


0.9130 


2 


10 


6 9609 


0.7980 




8 


4.0919 


0.9543 




10 


6.9615 


0.7980 




8 


4.1320 


0.9466 


2 


10 


7.1355 


0.7768 


■^2 


8 


4.3311 


0.8487 


^ 2 


10 


7.2010 


0.7224 




8 


4.3783 


0.8004 


^ 2 


10 


7.2489 


0.7168 




8 


4.3866 


0.7810 




10 


7.2543 


0.7631 


'Si 


8 


4.5628 


0.7656 


^F2 


10 


7.3103 


0.7568 


'Di 


8 


4.5642 


0.7654 


'P2 


10 


7.3563 


0.7046 


'G2 


8 


4.6678 


0.7885 


'Pi 


10 


7.4043 


0.6993 




8 


4.7087 


0.7220 


'Pi 


10 


7.5375 


0.6783 


'So 


8 


4.7653 


0.6836 


'F2 


10 


7.6303 


0.6757 


'So 


8 


4.8182 


0.6751 


'Fi 


10 


7.6670 


0.6720 


'Do 


8 


5.0062 


0.6472 


2 p 

"2,2 


10 


7.9283 


0.6407 


'S2,2 


8 


5.2186 


0.6140 


^^2,2 


10 


8.0142 


0.6331 



Figure captions 



Figure Al. Mapping in the ^A-parameter plane. 



Figure A2. The first excited triplet state for two electrons: 1, Energy versus (j) plot; 
2, the exact result; 3, energy value at the inflection point. 
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Figure A3. Spectrum for two electrons. 



Figure A4. Spectrum for three electrons. 



Figure A5. Spectrum for four electrons. 



Figure A6. Spectrum for five electrons. 



Figure A7. Ground state energy of quantum dot in magnetic field: (a) - 2 electrons, 
(b) - 3 electrons, (c) - 4 electrons, (d) - 5 electrons. Single electron ground state 
energies E{1) are subtracted. 




Fig. 1 (E. Anisimovas and A. Matulis "Energy Spectra 



3.9 





Fig. 3 (E. Anisimovas and A. Matulis "Energy Spectra ... ") 




Fig. 5 (E. Anisimovas and A. Matulis "Energy Spectra ... ") 




Fig. 6 (E. Anisimovas and A. Matulis "Energy Spectra ... ") 




Magnetic Field, Tesia 



Fig. 7a (E. Anisimovas and A. Matulis "Energy Spectra 
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Magnetic Field, Tesia 

Fig. 7b-d (E. Anisimovas and A. Matulis "Energy Spectra ... ") 



